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Abstract: We present results for the production of a Z boson in association with single top at 
next-to- leading order (NLO), including the decay of the top quark and the Z boson. This electroweak 
process gives rise to the trilepton signature + jets + missing energy. We present results for 

this signature and show that the rate is competitive with the contribution of the mixed strong and 
electroweak production process, tiZ. As such it should be observable in the full data sample from LHC 
running at y/s = 8 TeV. The single top + Z process is a hitherto unconsidered irreducible background 
in searches for flavour changing neutral current decays of the top quark in tt production. For a selection 
of cuts used at the LHC involving a &-tag it is the dominant background. In an appendix we also 
briefly discuss the impact of NLO corrections on the related tH process. 
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1 Introduction 

After only one year of 8 TeV running, the LHC has already become a tool for detailed studies of the 
top quark. With an increase to a higher centre-of-mass energy and anticipated integrated luminosities 
of up to 3000 fb _1 , the LHC will be able to achieve measurements of unprecedented precision in the 
top sector. With the advent of high statistics top physics, it will be possible to study not only the 
production of top quark pairs but also processes in which a vector boson is produced in association 
with top quarks. 

The CMS and ATLAS collaborations have produced first results on ttZ and ttW production in 
recent publications [1, 2]. The tiW process does not depend on the details of the top sector since the 
accompanying W boson is radiated from the initial state quarks. In contrast, the tiZ process directly 
probes the coupling of the Z boson to the top quark. Theoretical predictions are available for these 
processes at the NLO parton level [3-5] and in NLO calculations matched to a parton shower [6, 7]. 

In this context it is also interesting to consider the process where an extra Z boson is radiated in 
t-channel single top production. This predominantly proceeds through the leading order processes, 

u + b^d + t + Z, d + b^ u + t + Z , (1.1) 

for the production of a top quark, with smaller contributions from strange- and charm-initiated reac- 
tions. Production of an anti-top quark proceeds through the charge conjugate processes, 



d + b^u + t + Z, u + b^d + t + Z, 



(1.2) 




with a smaller rate at the LHC due to the difference in up- and down-quark parton distribution 
functions (pdfs). The leading order (LO) Feynman diagrams for the first process in Eq. (1.1) are 
shown in Fig. 1, including also the non- resonant contribution, diagram (g), that should be included 
when considering the charged lepton final state. The Z boson can be radiated from any of the four 
quark lines, or from the W boson exchanged in the i-channel. As can be seen from the diagrams, this 
process is related to hadronic WZ production by crossing. As a matter of principle, measurement of 
single top+Z is thus as important as measuring the W Z pair cross section, with the added bonus that 
it depends on the coupling of the top quark to the Z . In this paper, we present results for the single 
top + Z process to next-to-leading order (NLO) in QCD 1 . 

Although the single top + Z process is an electroweak one, in contrast to the QCD-induced pair 
production mode (ttZ), it contains fewer particles in the final state and is therefore easier to produce. 
Fig. 2 shows that any advantage in rate for the top pair production is effectively removed once an 
additional Z-boson is required. As a result, the single top + Z cross section is about the same size as 
the ttZ one. Given the status of current LHC searches for ttV production it is interesting to consider 
the expected experimental sensitivity to the single top + Z channel. In particular, the impact of these 
SM processes should already be present in current trilepton searches, albeit in regions of lower jet 
multiplicity. 

In order to properly assess the expected event rates in trilepton searches, in this paper we will 

1 Next-to- leading order QCD corrections to tZ associated production via the flavor-changing neutral-current couplings 
at hadron colliders have been considered in Ref. [8]. 
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Figure 2. NLO inclusive cross sections for single and top quark pair production with and without an accom- 
panying Z boson. The NLO ttZ cross section is estimated from the lowest order result using a if-factor of 
1.39 and renormalization and factorization scales fj, = mt + mz/2 [4]. 



consider the full process (and similarly for the charge conjugate process), 

u + b -> t + Z + d 

U fjT + fi+ (1.3) 
->■ v + e+ + b 

where the leptonic decay of the top quark is included and we have specified the charged leptons that 
are associated with the Z decay. The top quark decay is included using the techniques described in 
Refs. [9-11] and retains all spin correlations at the expense of requiring the top quark to be treated 
exactly on-shcll. Since this calculation involves an incoming 6-quark it is necessarily a five-flavor 
calculation. 

We have also considered the closely-related single top + H process which is of smaller phenomeno- 
logical interest in the Standard Model. A brief description of the next-to-leading order result is given 
in Appendix B. 
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2 Outline of calculation 



2.1 Leading order 

The leading order diagrams for this process are shown in Fig. 1. It is useful to consider the contribution 
from combinations of individual diagrams as follows: the Zj^f attached to the light quark line, M^ a ' b ', 
the Z/j* attached to the heavy quark line, M^- C ' d \ the Z/'y* attached to the t-channel W boson M^ e '^\ 
the non-resonant contribution with the lepton line attached to the i-channel exchanged W bosons, 
jlf(ff) . The computation of the amplitude can be performed in the unitary gauge. However a more 
compact expression is obtained in the Feynman gauge after the inclusion of an additional contribution 
representing the propagation of unphysical Higgs fields (represented by ip in diagram (f)). In the latter 
approach the cancellation of the terms associated with the longitudinal degrees of freedom is built-in. 
The explicit form of the leading order amplitudes is given in Appendix A. 

2.2 Next-to-leading order 

Next-to-leading order corrections to the single top + Z process are computed in a fairly straightforward 
manner. Virtual corrections to diagrams in which the Z boson is radiated from the i-channcl W or in 
which the lepton pair are produced in a non- resonant manner (c.f. Fig. l(e,f) and (g)) consist solely of 
vertex corrections and are therefore easily computed analytically. For the remaining diagrams, where 
the Z boson is radiated from one of the fermion lines, some of the vertex corrections can be computed 
in a similar fashion. However, the virtual amplitude also receives contributions from box diagrams 
containing three powers of the loop momentum. These corrections are computed numerically using a 
variant of the van Oldenborgh-Vermaseren scheme for the calculation of tensor integrals [12]. Scalar 
integrals are computed using the QCDLoop library [13]. We have also implemented a version of the 
usual Passarino-Veltman reduction algorithm [14], supplemented by special handling of regions of small 
Gram or Cayley determinants according to the procedure outlined in Rcf. [15]. In our implementation 
we find that the alternate reduction methods are used to improve the numerical stability of the 
calculation in approximately 0.3% of all events. 

As a further check, we compare the numerical calculation of the singular contributions to the 
amplitude to the known analytic form (after renormalization) [16], 
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where the invariants S25 and s\q are taken from the momentum assignment in equation (A. 7). The 
overall factor cr is, 

1 r(i + e)r 2 (i- e ) 

Cr= (47r)2-e r(l-2e) ■ (2 - 2) 

We find that less than 0.02% of all events fail this consistency check and are discarded. Moreover, 
these points lie in extreme phase space regions that contribute little to total cross sections. When 
realistic experimental cuts are applied the proportion of numerically unstable points removed from 
the calculation drops by a factor of about four. 

The calculation is performed in the four-dimensional helicity (FDH) scheme [17]. The mass renor- 
malization is fixed by the condition that the inverse propagator vanish on-shell. In the FDH scheme 
we have, 
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mw 


80.398 GeV 


IV 


2.1054 GeV 


m z 


91.1876 GeV 


Tz 


2.4952 GeV 


m t 


173.2 GeV 


G F 


1.116639 x 10~ 5 


a)f{m z ) 


0.130 


(4 h °{m z ) 


0.118 



Table 1. Input parameters used for the phenomenological results. The two values of as(mz) correspond to 
the choices made in the CTEQ6L1 and CTEQ6M pdf sets, used at LO and NLO respectively. 



and the wave function renormalization is, 



Zq = 1 — g 2 c r C F 



-+31n At 



+ 



(2.4) 



The coupling of the scalar ip to the quark field, proportional to the top mass, must also be renormalized 
in the same way. 

The top quark decay is included using the method of Rcf. [11]. We have included only the leading 
order amplitude for the decay since the rate for this process is already very small. 



3 Results 

For the results that we present in this paper, we have used the parameters listed in Table 1. From 
these, the Weinberg angle is fixed by the tree-level relation, 

S in 2 ^ = l-^, (3.1) 
m z 

which ensures that the amplitudes are gauge invariant. Since our calculation is performed in the 
five-flavour scheme, with an initial state massless 6-quark, we also set mi, = in the decay of the top 
quark. For simplicity we work in the framework of a unit CKM matrix. The parton distributions 
employed are the CTEQ6L1 set (used at LO) and CTEQ6M set (used at NLO) taken from rcf. [18]. 
The renormalization and factorization scales, denoted by /ir and [ip respectively, are taken to be the 
same for our standard scale choice, (1r — /if = rut- 

With these parameters, the total cross sections for tZ and tZ production as a function of the 
LHC operating energy ^fs are shown in Figure 3. Although the leading order process contains a 
quark, the i-channel exchange of the W boson means that the amplitude does not contain a collinear 
singularity and thus that the inclusive cross section is well-defined. The cross section for tZ production 
is approximately half the corresponding tZ rate, a reflection of the corresponding parton distribution 
function ratio, fd{x)/f u {x) ~ 0.5 at values of x typical of those relevant for this process, x > (m t + 
mz)/y/s ~ 0.02 — 0.03. The NLO corrections take a similar form for both processes, resulting in an 
increase in the cross section predictions of the order of 10%. Finally, we see that although the cross 
sections are only of the order of a few hundred femtobarns at yfs = 8 TeV, these processes have a 
combined cross section that is approximately a picobarn at y/s = 14 TeV. 

To investigate the scale dependence of this process we focus on the centre-of-mass energy -^/s = 
8 TeV. Since the tree level process does not contain a strong coupling the resulting cross section only 
depends on the factorization scale, but at next-to-leading order the renormalization scale enters for 
the first time. We find that varying both scales together in the same direction leads to an accidental 
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Figure 3. Cross sections for tZ and tZ production as a function of \fs. The leading order predictions are 
shown as dashed lines and the next-to-leading order solid lines. 



cancellation and therefore an artificially small estimate of the scale dependence. We therefore choose 
to vary them in opposite directions, = rm t , Hf = ni t /r with re [1/4,4]. The results are shown 
in Figure 4, where one can see that the overall scale dependence is still very weak. Even over such a 
large scale range the largest deviation from the central value is less than six percent. 

Before turning to less inclusive cases, we summarize our findings by presenting predictions for 
LO and NLO cross sections at y/s = 8 TeV. For the NLO prediction it is useful to consider the 
theoretical uncertainty that should be attributed to the calculation. In addition to the scale dependence 
uncertainty, based on the variation of r over the full range as described above, we also consider the 
effect of uncertainties in the extraction of the pdfs. By using the additional uncertainty sets provided 
in the CTEQ6 distribution, we find that this uncertainty is at the level of 7%. We thus find, 

<rLo(tZ) = 148 fb , cr NLO (tZ) = 160±2 (scale)+^ (pdf) fb , (3.2) 
a LO (tZ) = 68 fb , <r NLO {iZ) = 76t? (scalc)^ (pdf) fb , (3.3) 

Combining the two sources of error, the single top + Z cross section is thus predicted with a total 
uncertainty of just over 10%. 

3.1 Comparison of rates for tZ, iZ and ttZ 

As discussed in the introduction, the cross-section for ttZ production is comparable to that for the 
sum of tZ and iZ production. Referring to equation (1.3), the signature for tZ production is three 
charged leptons, missing energy (which can be reconstructed up to the usual two-fold ambiguity) and 
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Figure 4. Scale dependence of single top + Z cross sections at 8 TeV. The renormalization and factorization 
scales are varied in opposite directions according to fiR = rmt, /if = mt/r. 



jets. One of the jets may be 6-taggcd, although we ignore that possibility in this section. In the 
top-pair production scenario, the subsequent semi-leptonic decay of one top and the hadronic decay 
of the other, together with the leptonic decay of the Z-boson, gives rise to the same signature of three 
charged leptons, missing energy and jets. If some of the jets go undetected, then the question arises 
as to whether it is possible to disentangle these two production processes. 

In order to answer this question, we calculate jet-binned cross-sections for four processes, 

(a) v e e+b)Z , (c) v e e + b)i(-> qqb)Z , 

(b) t(-> e~v e b)Z , (d) t(-> qqb)t(-> e-y e b)Z , (3.4) 

with the decay Z — > fi~ n + understood in each case. We perform our comparison at the y/s = 14 TeV 
LHC. The scale \i = m t is used for the tZ and tZ calculations, and fx = m t + mz/2 for tiZ, following 
refs. [3, 4]. We will make use of three sets of kinematic cuts. The first, which we refer to as "standard 
cuts" , requires that the momenta of the leptons, jets and missing energy are each greater than 20 GeV, 
and that the pseudorapidity of the leptons and jets are constrained by \r)i\ < 2.5 and \r)j\ < 3.5. Jets 
are constructed with the anti-k t algorithm using AR = 0.4. The second set of cuts require a more 
central jet, \rjj\ < 2.0, but are otherwise the same. We shall refer to these cuts as "|?7j| < 2.0" cuts. 
The third set of cuts is identical to the standard cuts, but the jets are constructed using AR = 0.7. 
This is referred to as the "AR = 0.7" setup. 

The comparisons are shown for the LO results in figure 5. The figures on the left are for processes 
(a) and (c), which result in a final state signature with two positively charged leptons; the right-hand 
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Figure 5. Comparison of jet-binned cross-sections calculated at LO at y/s = 14 TeV. The left-hand plots 
show tZ production and tiZ production with the subsequent semi-leptonic decay of the top, resulting in a final 
state of fi~ fi + e + . The right-hand plots show iZ and tiZ production with the subsequent decay of the t, with 
a final state of [f /j, + e~ . The first row corresponds to the standard cuts described in the text, the second row 
uses the \rjj\ < 2.0 cuts, and the final row has the AR — 0.7 setup. The scale \i = mt is used for tZ and iZ, 
and = mt + mz/2 is used for ttZ. 



figures show processes (b) and (d), for which the signature includes two negatively charged leptons. 
Of course, the results for the tiZ process are the same irrespective of which top decays hadronically, 
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Jet multiplicity 





1 


2 


3 


Standard cuts 


LO 
NLO 


0.014 
0.011 


0.331 
0.237 


1.05 

0.585 


0.693 


tel < 2 


LO 
NLO 


0.140 
0.115 


0.856 
0.669 


0.400 
0.531 


0.211 


AR = 0.7 


LO 
NLO 


0.014 
0.010 


0.336 
0.241 


1.05 
0.661 


0.614 



Table 2. Jet-binned LO and NLO cross-sections (in fb) for tZ + tZ production at the ^/s = 14 TeV LHC, for 
the three sets of cuts described in the text. 

whereas the tZ cross-sections are a factor of approximately two greater than those for iZ, as indicated 
in figure 3. This feature suggests a way of distinguishing between the single top + Z and tiZ processes, 
by, for example, considering the asymmetry between l + l~l' + and production. This method 

would rely on a stringent rejection of backgrounds, some of which would display a similar asymmetry. 

The first row of figure 5 corresponds to the standard set of cuts. It is seen that most of the jets 
in tZ production arc able to pass these cuts, so that the two-jet bin dominates the total cross-section. 
By contrast, the tiZ process has a small cross-section in the two-jet bin and a negligible contribution 
to the one-jet bin. 

The effect of lowering the cut on the jet pseudorapidity to \rjj\ < 2.0 is shown in the second 
row. Since one of the jets in tZ production is usually quite forward, with the other one central, it 
is unsurprising to see that the one-jet bin is dominant for tZ production. It is also evident that the 
stricter jet cut has shifted some of the ttZ events to the lower jet bins, with the result that the two-jet 
bin contains a significant proportion of events originating from this process. 

The third row shows the results using the AR. = 0.7 setup. This has little effect on the jets 
originating from tZ production: since one is forward and the other one central, there is little oppor- 
tunity for these to be clustered into one jet. The effect is more pronounced for ttZ , enhancing the 
cross-section in the two-jet bin. 

The effect of NLO corrections to the tZ + iZ cross-sections are shown in table 2. The total cross- 
section shows a slight increase from a^o = L4 fb at LO to <jnlo = 1-5 fb at NLO. However, looking 
at the standard cuts, it is clear that this increase is not uniform over the jet bins. The three-jet bin 
contributes around half the total cross-section, indicating that the additional radiated gluon is usually 
quite hard. This has the effect of migrating events from one jet bin to the next, with the result that the 
cross-sections in the zero-, one- and two-jet bins decrease due to the NLO corrections. This holds true 
when a larger jet is used, AR = 0.7, although the two-jet bin is larger and the three-jet bin smaller 
than with the standard cuts. This is because of the increased likelihood of clustering the radiated 
gluon with one of the LO partons, leaving two jets. When the \r]j\ < 2.0 cuts are used, the NLO 
corrections decrease the one-jet bin and increase the two-jet bin. The three-jet bin is much smaller 
than for the standard set of cuts. 

The NLO corrections indicate that distinguishing between tZ and tiZ production may be more 
difficult than a LO calculation leads one to expect. The NLO corrections deplete the tZ cross-sections 
in the bins where they are dominant over the tiZ cross-sections, and result in comparable cross- 
sections in the three-jet bin, which only received contributions from ttZ at LO. Nor is this the final 
story. A more realistic calculation of the jet-binned cross-sections would take parton showering into 
account. These effects can have a significant impact on exclusive observables. It should also be borne 
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in mind that NLO corrections and/or parton showering effects may modify the tiZ results. Ideally, a 
comparison would be performed after calculating both processes to NLO in QCD, and then interfacing 
them with a parton showering program that preserves the NLO accuracy. 

3.2 Single top + Z as a background in non-standard top decay searches 

The top quark decays primarily via a W boson, t — > Wq, with a bottom quark being the most 
likely decay product and the presence of strange or down quarks suppressed by the off-diagonal CKM 
elements. In the standard model, decays through a flavor-changing neutral current (FCNC) are loop- 
suppressed, yielding a very small branching ratio B(t — > Zq) < 10~ 12 [19]. Therefore, the observation 
of such a decay would be indicative of New Physics. Searches for FCNC decays in ti production were 
conducted by both CDF [20] and DO [21]. Currently, the best constraints come from ti production at 
the LHC: ATLAS constrains the branching ratio B(t -> Zq) < 0.73% with 2.1 ftr 1 of data at y/s = 7 
TcV [22], while CMS constrains B(t -> Zq) < 0.24% with 5.0 ftr 1 of data at the same energy [23]. 

As the second top is taken to decay through the Standard Model mode t — > Wb, the signature 
of these events (with leptonic decays of both the W- and Z-bosons) is three charged leptons, missing 
energy from a neutrino (whose longitudinal momentum is reconstructible, up to the usual two-fold 
ambiguity), and two or more jets, one of which can be 6-tagged. The same signature is expected in tZ 
and tZ production. However, neither the ATLAS [22] nor the CMS [23] analysis take this background 
into account. The purpose of this section is to look at the role of tZ and iZ production as a background 
to FCNC top decays. 

We consider decays of the W- and Z bosons into different flavored leptons, Z —> /J~/i + and 
W — > v e e : and impose a set of cuts similar to those used in the CMS analysis 2 : 

• Leptons are required to have transverse momentum prj > 20 GeV and pseudorapidity |rj;| < 2.5. 

• The missing transverse momentum is constrained by pr.miss > 30 GeV. 

• Jets are defined with the anti-fey algorithm with AR = 0.5, and are required to have pT.j > 30 
GeV and \rjj \ < 2.4, and to be separated from any lepton by ARji > 0.4. 

• The same- flavor dilepton pair is required to have mass 60 GeV < mu < 120 GeV. This pair is 
taken as originating from the Z-boson, with the remaining lepton originating from the VF-boson. 

• Each lepton is required to be isolated. In particular, the ratio of the sum of the transverse 
energies and momenta of all objects (leptons and jets) within AR = 0.3 of the lepton to the 
lepton's transverse momentum must be less than 0.125 for leptons originating from the Z-boson, 
and less than 0.1 for the lepton originating from the VF-boson: 

Ear w <o. 3 ( e t+Pt) Y,ar z <o. 3 ( e t+Pt) 
Ptj Pt.i 

(for our purposes, we set Ex = Px)- 

In addition to the above cuts, CMS uses two further sets of cuts, called "St" cuts and "6-tag" 
cuts. In the case of the former, the following cuts are applied: 

• At least two jets are required, with the transverse momentum cut as above. 

• The total transverse momentum Sx = Ylj PT.j + Yli Pt.i + PT.miss > 250 GeV. 

2 The cuts used by CMS are slightly more complicated, since they take into account various detector effects. 
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St cuts 


6-tag cuts 


Ztj 


onlo 


33.3(l)±|;g 
52.0(1)^;| 


i4.3(i)^: 6 8 
24.5(i);?.| 


Ztj 


&NLO 


17.5(l)t? : 6 o 
26.2(1);°:! 


7.7i(i)±° : i| 

12.5(l)+g;| 



Table 3. Leading- and next-to-leading order cross-sections (in ab) for Z(—t /i~/i + )i(— > v e eb)j using the two 
classes of cuts used in the CMS searches for FCNC in top decays. The cross-sections are evaluated at a 
scale /i = mt, with the integration error in the last digit in parentheses. The effect of using a scale choice of 
fi = mt/2 and fi = 2mt are shown as subscripts and superscripts respectively. 

• The masses of the Zj and 1^6-system are constrained to be between 100 GeV and 250 GcV. 
The "6-tag" cuts are: 

• At least two jets are required, one of which is 6-taggcd. 

• The masses of the Zj- and W6-systcms are constrained to be close to the top mass: \mzj — 7Tit| < 
25 GeV and \mwb — mt\ < 35 GeV. 

The LO and NLO cross-sections for tZ and iZ production are shown at the y/s = 7 TcV LHC in 
table 3. There is a negligible change when the three charged leptons have the same flavor. We note 
that the NLO corrections have a substantial effect on the cross-sections, with a K- factor of around 1.5 
for the St cuts and 1.7 when the ^-tagging cuts are used. This is because the additional jet from the 
real radiation helps to satisfy the jet cuts. The scale uncertainty is larger than discussed previously 
for the inclusive production, and we estimate these uncertainties by varying both the factorization and 
rcnormalization scales in the same direction, between mt/2 and 2mt- This gives a scale uncertainty 
of around 5-7%. The pdf uncertainty is not taken into account, but is expected to be similar in 
magnitude. 

The dominant background in the CMS analysis comes from WZjj production, with leptonic decay 
of the weak bosons. Imposing the St cuts we calculate this cross-section to be 0.91 fb at LO, with 
a scale uncertainty of around 25%. Multiplying by a factor of four to include all leptonic final states 
eee, ee/i, /x/ie, n^n, we find in a sample of 5.0fb~ 1 that this corresponds to 0.91 x4x5 = 18.2 events. 
This is consistent with the CMS calculation of 13.6 ± 2.6 WZjj events. 

We can convert the cross sections of table 3 into event rates to compare with the CMS study in 
similar fashion. This implies that 1.6 events should be seen for the tZ + iZ background when the St 
cuts are used. This is a small but not negligible increase on the 16.2 overall background events that 
are expected. However, when the 6-tag cuts are used, the overall CMS background estimation drops 
significantly to 0.83 events, due to a more stringent cut on the mass-window of the weak boson-jet 
system, and the requirement of a 6-tag. Since our implementation of tZ production constrains the 
Wfr-system to the top mass and guarantees the presence of a 6-jct, the effect of these cuts is far less 
severe, and we expect 0.74 events coming from the tZ + tZ background with this set of cuts. At 
present, the best constraint on the FCNC branching ratio is found using the St cuts. However, it is 
possible that this situation could be changed once the dominant single top + Z contribution to the 
backgrounds with 6-tag cuts is included. 
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4 Conclusions 

We have calculated the production cross-section of single top + Z-boson to NLO in QCD, including 
the leptonic decays of the top quark. We have demonstrated that this process is competitive in rate 
with the mixed strong and electroweak ttZ process. As such, it should be observable in recorded 
data from the LHC, despite being subject to a considerable reducible background from W ± Z+2 jet 
processes. Given this, the potential to constrain the top-Z boson coupling through the tZ process 
should be investigated further. Moreover, we have shown that the use of jet-binned cross-sections may 
be helpful in distinguishing this process from the ttZ process, although this requires further effort 
on the theoretical front to determine the effects of parton showering for this observable. In addition, 
this process constitutes an irreducible and potentially dominant background in searches for flavour 
changing neutral current decays in ti production, which is not taken into account in current searches. 
It will be challenging to remove because, like the signal, it contains a real top quark. Code for this 
phenomenological interesting process, as well as the related tH process, is included in MCFM v6.6. 
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A Calculational details 

A.l Notation for spinor products 

We adopt the following notation for massless spinors, 

\i) = \i+) = u+{pi), \i] = \i-) = u-(pi) , 

(*| =(*-! = u-(Pi)> [*l = (* + I = S +CP*) ■ (A.l) 
Further the spinor products are defined as, 

(ij) = (i ~ \j+) = u^(p l )u+{p :j ) , 

[ij] = (i+ = u + { Pi )u-( Pj ) , (A.2) 
with pi, pj massless particles. With our convention, 

(ij) \ji] = 2 Pi -P] = s ij ■ ( A - 3 ) 

We shall use the standard trick [24] of decomposing the massive momentum, p 2 = m 2 into the 
sum of two massless momenta, p = p^ + ar\ with the constant a given by, 



TO? 



(vMv] ' 

We may write the massive spinors as combinations of massless spinors as follows, 

u-(p) = M(y+ m t)i^TT> ft+(p) = (Vt\(jf + m) , b , . (A.5) 
[VP I (VtP ) 

v+(p) = (^-m t )[v)y—r, u-(p) = (jf - mt)\v]y^j\ ' ( A ' 6 ) 

The spin labels of the massless spinors encode the polarization information of the massive 

quarks and they are equivalent to helicities only in the massless limit. 
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A. 2 Lowest order matrix element 

We present results for the basic amplitude at leading order, 

u(pi) + b(p 2 ) -t l{jpa) + a(pi) + t(p 5 ) + d(p 6 ) , 



(A.7) 



where I, a are the lepton and anti-lcpton respectively and momentum labels for the particles are given 
in parentheses. 

We begin by introducing the relevant couplings that appear in the calculation. The current for 
the emission of a Z boson or virtual photon that decays into a left-handed lepton pair enters with a 
strength, 

Vj = QjQe + LjhsuDzissA) , Vf = Qjq e + RjleSsiDzisu) , (A.8) 

where the superscript denotes the helicity of the outgoing quark and the subscript the flavor of the 
quark from which the boson is emitted (j = u,d). In this formula the individual quark and lepton 
couplings are themselves defined by, 



h = 



le 



tj — 2Qj sin 2 9w 

sin 28w 
- 1 — 2q e sin 2 9\y 



Rj = 



—2Qj sin 2 9w 

sin 26w 
-2q e sin 2 9 W 



where q e 



sin 29w sin 29y/ 

T, t u = 1 and = — 1. The Z propagator denominator is, 

Dz{s3i) = o" ■ 

S34 - m z 



(A.9) 
(A.10) 

(A.ll) 



We first consider the case of a negative helicity outgoing lepton and a negative spin-label for the 
top quark. The contributions to the amplitudes, calculated in the Feynman gauge and labelled by the 
diagrams in Fig. 1 are, 



(1- , 2 6 - , 3," , 4+ 5 t " , 6+) = D w (s 25 ) — 

(5 b 6\ [1 4] (311 + 412] - ^ (36) [1 2] (5 b |3 + 6|4 

Sl34 \ / S 346 \ 



(A.12) 



MM)(l-2^3 r ,4+5 t -6+) = ^M 



-^(3|(2 + 4)|l](65 b \[24] 
S234 \ / 



S34 
V L 



M^)(l-,2 b -,3r,4+ 5 t " 6+) = 



(s 345 - m 2 ) 
Dw(s 2 5)D w (s w ) 



V«ml (3 6) [12] [4 V ] 

(s 3 45 - TO 2 ) [5 b T]) 

(6|(1 + 2)|4] (35 b \[12] 



(A.13) 



•S34 



-(Vu-v d L ) 



{ (3|(1 + 6)|4] (6 5 b ) [12] + ^5 b |(1 + 6)|2] (3 6) [14] + (6|(3 + 4)|1] ^3 5 b ) [2 4] } 

j^M{«?-Kf-v- + l ?}" 



m|_ (3 6) [14] [2 7?] 
2 



M^(l-, 2^,3^,4+5,-, 6+) 



Dw{s25)Dw(sw) 
2 sin 2 9 W S235 



3 5 b )[14](6|(l + 4)|2] 



(A.14) 
(A.15) 
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For the case of a positive spin-label for the top quark we have, 



M<»' 6 >(1- 2 6 " 3j-,4+ 5+ 6+) = D w (s 25 )^ 

S34 



V U L (6 t,) [14] 



(3|(1+4)|2] 



(3 6) [12] 



S134 (5 b rj) s 346 (5 b 77) 

M^\ K ^, K Ai^^)- Dw{sw)mt 

V d L <3|(2 + 4)|1] <677>[2 4] 



(r?|(3 + 6)|4] 



S34 



(3 6) [12] 



(S345 - "i?) 
(6|(1 + 2)|4] (377) [12 



4 5 c 



S234 



(S345 - mf) 



(A.16) 



(A.17) 



M^)(l-,2 6 -,3 ; - 4+ 5+ 6+) = D w (s 25 )D w ( Sl6 ) 



S34 



Vj 



{ 1 «3|(1 + 6)|4] < 6? ?> [1 2] + <77|(1 + 6)|2] (36) [14] + (6|(3 + 4)|1] (3ry) [24]} 



(5 b 77) 
5 (3 6) [14] 



2 5 c 



[Vu - v d L + V U H - V d R } 



M^)(l-,2 6 -,3r,4+ 5+6+) 



^(S25)^(si 6 )m t r(3 77) [14] (6|(1 + 4)|2] 



2 Sin 6»iyS 23 5 



<5 b v) 



(A.18) 
(A.19) 



Note that the opposite helicity combination for the lepton line is obtained by performing the flip 
3H4,l e 4r e for M ( a ' b \ M << c > d \ M^ eJ \ The amplitude docs not contribute for the opposite 

helicity. 

The total leading order amplitude is obtained by summing these four subamplitudcs. In order 
to allow the Z boson to be off-shell but still retain gauge invariance, we use a simple prescription to 
incorporate the Z width [25]. We use the propagator factor Dzis-a) in the amplitudes as written 
above and then multiply the whole amplitude by, 



S34 - rn z 



s 34 - to| + imzTz 



(A.20) 



B Associated production of a single top and Higgs boson 

In this appendix we briefly describe the NLO calculation of single top + Higgs boson production, which 
is very similar in many respects to the single top + Z process that is the main topic of this paper. In 
the limit in which the light quarks are taken to be massless, there are only two leading order diagrams, 
as shown in Figure 6, with the Higgs boson attaching to either the top quark or the t-channel W boson. 
This process has previously been considered in Refs. [26-28]. The gauge cancellation between the two 
diagrams in Fig. 6 results in a smaller cross section compared to the associated pair production mode, 
ttH. In addition, because of the small branching ratios of a 126 GeV Higgs boson to the cleanest 
modes (H — > four leptons and H — !> 77), single top + H production will be extremely challenging to 
observe. Nevertheless, like the ttH process, this channel has the potential to measure the coupling 
of the Higgs boson to the top quark. Reliable theoretical estimates for the ttH process, accurate to 
NLO, are given in Refs. [29-32] and including also the effect of a parton shower in Refs. [33, 34]. Here 
we bring the accuracy of the single top + H channel to the NLO parton level. 

Our results are calculated using the same numerical procedure described in Section 2. Due to 
the simplicity of the scalar coupling of the Higgs, it is possible to immediately reduce the rank of the 
tensor integrals that appear in the 1-loop calculation to a maximum of two. As a result we find that 
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(a) (b) 

Figure 6. Feynman graphs to calculate the lowest order amplitudes for single top + Higgs production. The 
diagrams where the Higgs boson couples to the light quarks have been dropped. 



the calculation is significantly more stable than the single top + Z case, with an order of magnitude 
less events discarded due to insufficient numerical precision in the pole terms (less than 0.005%). The 
renormalization of the Yukawa coupling of the Higgs boson to the top quark takes exactly the same 
form as the renormalization of the tp coupling already discussed in Section 2. 

For the results presented here we use mjy = 126 GeV based on the first observation of a new 
boson at the LHC. The cross sections for tH and tH production as a function of the LHC operating 
energy yfs are shown in Figure 7 (left). The effect of next-to- leading order corrections is larger than 
in the single top + Z case, with an increase in the cross section of approximately 15% at NLO. 
To investigate the scale dependence of this process we focus on the case yfs = 8 TeV. In contrast 




Vs [TeV] ^ m t 

Figure 7. Left: Cross sections for tH and iH production as a function of y/s. Right: Scale dependence 
of single top + H cross sections at 8 TeV (fj. = fiR = hf)- In both cases, the Higgs boson has a mass 
uih = 126 GeV and leading order predictions are shown as dashed lines, next-to-leading order as solid lines. 

to the production of single top + Z, in this case we find the largest scale dependence when both 
renormalization and factorization scales arc varied together. The results are shown in Figure 7 (right), 
where we consider scale variation by a factor of four about the central value, \x = m t . Once again the 
NLO scale dependence is very mild, as expected in an electroweak process. 

This process has received considerable interest recently as a probe of non-standard couplings of 
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the Higgs boson to top quarks [28, 35]. If the couplings deviate from their SM values (e.g. due to 
New Physics effects in loops) then the tH cross-section may be significantly enhanced. We allow the 
possibility of anomalous couplings in our code to enable a NLO calculation of such effects. 



References 

[1] CMS Collaboration, Measurement of the Associated Production of Vector Bosons with Top-Antitop 
Pairs at 7 TeV, Tech. Rep. CMS-PAS-TOP-12-014, CERN, July, 2012. 

[2] Search for tiz production in the three lepton final state with 4.7 fb of i/s = 7 tev pp collision data 
collected by the atlas detector, Tech. Rep. ATLAS-CONF-2012-126, CERN, Geneva, Aug, 2012. 

[3] A. Lazopoulos, T. McElmurry, K. Melnikov and F. Petriello, Next-to-leading order QCD corrections to 
tiZ production at the LHC, Phys.Lett. B666 (2008) 62-65 [0804.2220]. 

[4] A. Kardos, Z. Trocsanyi and C. Papadopoulos, Top quark pair production in association with a Z-boson 
at NLO accuracy, Phys.Rev. D85 (2012) 054015 [1111.0610]. 

[5] J. M. Campbell and R. K. Ellis, tiW ± production and decay at NLO, J HEP 1207 (2012) 052 
[1204.5678]. 

[6] M. Garzelli, A. Kardos, C. Papadopoulos and Z. Trocsanyi, Z0 - boson production in association with a 
top anti-top pair at NLO accuracy with parton shower effects, Phys.Rev. D85 (2012) 074022 
[1111.1444]. 

[7] M. Garzelli, A. Kardos, C. Papadopoulos and Z. Trocsanyi, t tbar W and t tbar Z Hadroproduction at 
NLO accuracy in QCD with Parton Shower and Hadronization effects, 1208.2665. 

[8] B. H. Li, Y. Zhang, C. S. Li, J. Gao and H. X. Zhu, Next-to-leading order QCD corrections to tZ 
associated production via the flavor-changing neutral- current couplings at hadron colliders, Phys.Rev. 
D83 (2011) 114049 [1103.5122]. 

[9] R. Kleiss and W. J. Stirling, Top quark production at hadron colliders: some useful formulae, Z.Phys. 
C40 (1988) 419-423. 

[10] J. M. Campbell, R. K. Ellis and F. Tramontano, Single top production and decay at next-to -leading 
order, Phys.Rev. D70 (2004) 094012 [hep-ph/0408158]. 

[11] J. M. Campbell and R. K. Ellis, Top-quark processes at NLO in production and decay, 1204.1513. 

[12] G. van Oldenborgh and J. Vermaseren, New Algorithms for One Loop Integrals, Z.Phys. C46 (1990) 
425-438. 

[13] R. K. Ellis and G. Zanderighi, Scalar one-loop integrals for QCD, JHEP 0802 (2008) 002 [0712.1851]. 

[14] G. Passarino and M. J. G. Veltman, One Loop Corrections for e+ e- Annihilation Into mu+ mu- in the 
Weinberg Model, Nucl. Phys. B160 (1979) 151. 

[15] A. Denner and S. Dittmaier, Reduction schemes for one-loop tensor integrals, Nucl. Phys. B734 (2006) 
62-115 [hep-ph/0509141]. 

[16] S. Catani, S. Dittmaier and Z. Trocsanyi, One loop singular behavior of QCD and SUSY QCD 
amplitudes with massive partons, Phys.Lett. B500 (2001) 149-160 [hep-ph/0011222]. 

[17] Z. Bern, A. De Freitas, L. J. Dixon and H. Wong, Supersymmetric regularization, two loop QCD 
amplitudes and coupling shifts, Phys.Rev. D66 (2002) 085002 [hep-ph/0202271]. 

[18] J. Pumplin, D. Stump, J. Huston, H. Lai, P. M. Nadolsky et. at, New generation of parton distributions 
with uncertainties from global QCD analysis, JHEP 0207 (2002) 012 [hep-ph/0201195]. 

[19] E. N. Glover, F. del Aguila, J. Aguilar-Saavedra, M. Beccaria, S. Bejar et. al, Top quark physics at 
colliders, Acta Phys.Polon. B35 (2004) 2671-2694 [hep-ph/0410110]. 

[20] CDF Collaboration, T. Aaltonen et. al., Search for the Flavor Changing Neutral Current Decay t — > Zq 
mpp Collisions at sfs = 1.96 TeV, Phys. Rev. Lett. 101 (2008) 192002 [0805.2109]. 



- 16 - 



[21] DO Collaboration, V. M. Abazov et. al, Search for flavor changing neutral currents in decays of top 
quarks, Phys.Lett. B701 (2011) 313-320 [1103.4574]. 

[22] ATLAS Collaboration, G. Aad et. al, A search for flavour changing neutral currents in top-quark 
decays in pp collision data collected with the ATLAS detector at sqrt(s) = 7 TeV, JHEP 1209 (2012) 
139 [1206.0257]. 

[23] CMS Collaboration, S. Chatrchyan et. al., Search for flavor changing neutral currents in top quark 
decays in pp collisions at 7 TeV, 1208.0957. 

[24] R. Kleiss and W. Stirling, Spinor Techniques for Calculating pp — > W ± /Z° + Jets, Nucl.Phys. B262 
(1985) 235-262. 

[25] U. Baur, J. Vermaseren and D. Zeppenfeld, Electroweak vector boson production in high-energy ep 
collisions, Nucl.Phys. B375 (1992) 3-44. 

[26] F. Maltoni, K. Paul, T. Stelzer and S. Willenbrock, Associated production of Higgs and single top at 
hadron colliders, Phys.Rev. D64 (2001) 094023 [hep-ph/0106293]. 

[27] V. Barger, M. McCaskey and G. Shaughnessy, Single top and Higgs associated production at the LHC, 
Phys.Rev. D81 (2010) 034020 [0911.1556]. 

[28] M. Farina, C. Grojean, F. Maltoni, E. Salvioni and A. Thamm, Lifting degeneracies in Higgs couplings 
using single top production in association with a Higgs boson, 1211.3736. 

[29] W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira et. al., Higgs radiation off top quarks at 
the Tevatron and the LHC, Phys.Rev. Lett. 87 (2001) 201805 [hep-ph/0107081]. 

[30] L. Reina, S. Dawson and D. Wackeroth, QCD corrections to associated t anti-t h production at the 
Tevatron, Phys.Rev. D65 (2002) 053017 [hep-ph/0109066]. 

[31] S. Dawson, C. Jackson, L. Orr, L. Reina and D. Wackeroth, Associated Higgs production with top quarks 
at the large hadron collider: NLO QCD corrections, Phys.Rev. D68 (2003) 034022 [hep-ph/0305087]. 

[32] W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira et. al, NLO QCD corrections to t anti-t 
H production in hadron collisions, Nucl.Phys. B653 (2003) 151-203 [hep-ph/0211352]. 

[33] M. Garzelli, A. Kardos, C. Papadopoulos and Z. Trocsanyi, Standard Model Higgs boson production in 
association with a top anti-top pair at NLO with parton showering, Europhys.Lett. 96 (2011) 11001 
[1108.0387]. 

[34] R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, R. Pittau et. al., Scalar and pseudoscalar Higgs 
production in association with a top-antitop pair, Phys.Lett. B701 (2011) 427-433 [1104.5613]. 

[35] S. Biswas, E. Gabrielli and B. Mele, Single top and Higgs associated production as a probe of the Hti 
coupling sign at the LHC, 1211.0499. 



- 17 - 



